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THEORETICAL STABILITY AND CONTROL CHARACTERISTICS OF WINGS WITH 

VARIOUS AMOUNTS OF TAPER AND TWIST 

By Henry A. Pearsow and Robert T. Jones 



SUMMARY 

Stability derivatives have been computed for twisted 
wings of different flan forms that include variations in 
both the wing taper and the aspect ratio. Taper ratios 
of 1.0, 030, and 0.25 are considered for each of three 
aspect ratios: 6, 10, and 16. The specific derivatices for 
which results are given are the rolling-moment and the 
yawing-moment derivatives with respect to (a) rolling 
velocity, (b) yawing velocity, and (c) angle of sideslip. 
These results are given in such a form that the effect of 
any initial symmetrical wing twist (such as may be 
produced by flaps) on the derivatives may easily be taken 
into account. 

In addition to the stability derivatives, results are 
included for determining the theoretical rolling moment 
due to aileron deflection and a series of influence lines is 
given by which the loading across the span may be de- 
termined for any angle-of-attack distribution that may 
occur on the wing plan forms considered. The report also 
includes incidental references to the application of the 
results. 

INTRODUCTION 

Although a formal theory for the dynamics of air- 
plane motions has been available for many years, air- 
plane designers have not been in a position to utilize 
this theory to its fullest advantage on account of lack 
of knowledge of the basic physical quantities involved. 
It is true that the physical quantities, or stability 
derivatives, have been determined by test or calcula- 
tion in a number of instances, but there exists no sys- 
tematic series or correlation of tests sufficient to guide 
the designer in the prediction of these factors. 

As is well known, the calculations involved in aero- 
dynamic wing theory have been developed and refined 
to such an extent that it is possible to predict quite accu- 
rately the air moments and forces on the isolated wing 
at a fixed speed and incidence. Since several of the 
airplane lateral-stability derivatives depend almost 
entirely on the aerodynamic characteristics of the wing 
and since it would be desirable in any case to know the 
separate effects of variation of wing form on stability, 
it was thought worth while to extend the calculations 
to the determination of the moments developed by the 
wings when the airplane is disturbed from steady flight. 



This report gives theoretical stability derivatives for 
a variety of wing shapes including nine different plan 
forms and covering, in most cases, an arbitrary distri- 
bution of twist. 

Past work on the stability characteristics of wings has, 
except in isolated cases, been confined to analysis by the 
"strip method," wherein the effects of aerodynamic in- 
duction were neglected. The main effects of the induc- 
tion are included in the present computations, although 
the secondary influence of distortion or curvature of the 
wake is neglected. 

DEFINITIONS 

The axes used in specifying moments, angular veloci- 
ties, etc., are fixed in the wing and therefore move rela- 
tively to the air and to the earth. The X axis passes 
through the wing aerodynamic center in the plane of 
symmetry and is so chosen as to point directly into the 
line of the relative wind when the wing is moving stead- 
ily. Otherwise the axes form an orthogonal system as 
shown in the back cover of the report. 

The derivatives that may be obtained enable an esti- 
mate to be made of the variation of both rolling moment 
and yawing moment with (1) rolling velocity, (2) yaw- 
ing velocity, and (3) sideslip angle. These factors, des- 
ignated by Ct v , C Ut , etc., are to be used in the following 
general formulas to determine the wing rolling and yaw- 
ing moment in combined rolling, yawing, and sideslip- 
ping motion: 

^=<fr)+°<lp)+^ « 

Subscripts p and r are used to designate the partial 
derivatives of the well-known wing rolling-moment and 
yawing-moment coefficients, C t and (?», with respect to- 
instantaneous rolling and yawing angular velocities (ex- 
pressed as helix angles) and /3 is used to designate the 
partial derivatives of these coefficients with respect to 
instantaneous sideslip angles. In this manner the no- 
tation is considerably shortened from the usual more 

cumbersome expressions bCifb (j^fr)' (^"p-)' etc - 

Expressing the rolling and yawing moments as the sums- 
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of partial linear factors is considered valid for motions 
that are slow relative to the flight speed V and for small 
displacements, such as occur in ordinary unstalled 
maneuvers and such as are considered in the study of 
stability. 

a, angle between the zero-lift direction of the wing 
section and the air velocity at infinity, radians. 

8, parameter defining spanwise position, y= — | cos 6 

(when 6=0, y=~;d=ir, 2/=|)- 

Co, C}, d, coefficients of cosine series expressing 

wing plan form. 
Oi, rolling-moment coefficient.- 
C a , yawing-moment coefficient. 
p, angular velocity in roll, radians per sec. 
r, angular velocity in yaw, radians per sec. 
V, flight velocity of wing along X, f . p. s. 
i8, angle of sideslip, radians. 
5, aileron deflection, radians. 
C, pi rate of change of roUing-momenfr^coeflicient C t 

with the helix angle pb/2V. 
C Xp , rate of change of yawing-moment coefficient C„ 

with the helix angle pb/2V. 
Ci T , rate of change of rolling-moment coefficient C t 

with the helix angle rb/2V. 
C» r , rate of change of yawing-moment coefficient C n 

with the helix angle rb/2V. 
Ci0, rate of change of rolling-moment coefficient C t 

with sideslip angle /S. 
C K p, rate of change of yawing-moment coefficient C„ 

with sideslip angle 0. 
C ts , rate of change of rolling-moment coefficient Ci 

with aileron angle 5. 
C„ a , rate of change of yawing-moment coefficient C, 
with aileron angle 5. 
L, total wing rolling moment, ft. -lb. 
N, total whig yawing moment, f t.-lb. 
q, dynamic pressure, lb. per sq. ft. 

5, wing area, sq. ft. 

c, chord length at any section, ft. 
c„ chord length at plane of symmetry, ft. 
m 0 , section slope of the lift curve, per radian. 

6, wing span, ft. 

A, wing aspect ratio, br/8. 

A n , B n , coefficients of Fourier series. (See 

reference 2.) 
C|, section lift coefficient, section lilbjqcdy. 
C L , wing lift coefficient, wing lift/gS'. 
c d0 , section profile-drag coefficient. 
c i( , section induced-drag coefficient. 
Odq, wing profile-drag coefficient. 
X, taper ratio: i. e., ratio of the fictitious tip chord, 
obtained by extending the wing leading and 
trailing edges to the tip, to the root chord, 
r, dihedral angle, radians. 



CONDITIONS RELATING TO THE COMPUTATIONS 

PLAN FORMS 

The particular chord distributions for which the com- 
putations were made are illustrated in figure 1. Table 
I gives the coefficients of the cosine series used to express 
these chord distributions in terms of 8. Although the 
quarter-chord line is shown to bo straight, it is permis- 
sible to apply the results to wings with similar chord 
distributions but with the different plan forms that 




(a) ,4-8. 

(b) .4=10. 
(0) A-10. 

Figure 1.— Wing plan (onus considered. 



may be obtained by small alterations of the shapo of 
the quarter-chord line. The computations were made 
for three aspect ratios, 6, 10, and 16, and for threo taper 
ratios, 1.00, 0.50, and 0.25. The wing plan forms used 
only approximate those of linearly tapered wings with 
rounded tips. 

LIFT DISTRIBUTIONS 

Kolling, yawing, and sideslipping motions introduce 
varying resolutions of the relative-wind velocity over 
the wing. It is evident that these variations can, to a 
certain extent, be replaced by a fictitious warp or twist 
of the wing in straight flight. The procedure followed 
here is to calculate the spanwise lift and drag distribu- 
tions for the fictitious twist (i. e., that replacing the 
effect of motion) in the ordinary way, but to incline the 
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The values of G R ^fa s given in figures ,12 and 13 

have "been found to be in error because of the omnission 
of a tern in the expansion of the formula for the yawing 
derivative due to yawing. The correct values for these 
figures, which have been supplied by LHAL , nay be obtained 
fron the following table: 



VALUES 02 C n /a s 



\. Extent of flap 
! from wing 
■^v. center 

A 


' 0.25 


0.50 


0.75 


1.00 


X = 0.25 


6 
10 
16 


a 0.O06 
a .008 
a .005 


0.058 
.052 
.039 


a 0.175 
a .135 
a .C98 


0.290 
.220 
.147 


X = 0.50 


6 f a 0.0O2 
10 'j a .CQ4 
16 | a .003 

i 


0.060 
.050 
.039 


a 0.208 
a .l65 
a .120 


0.364 
.284 
.211 


X = 1.00 


10 
■ IS 


a -0.0Cl 
a .001 
a -.001 


•G.059 
.044 
. .035 


a 0.240 
a .182 
a .149 


0.450 
.365 
.297 


Elliptical wing 


10 , | 0.006 I 0.050 ! 0.142 


" 0.242 



Obtained fron faired curves. 
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Relative distance from wing center 

(a) X=L00. Cb] W1JS0. Co) X-0.25. 

Ficcbi 2.— Load distribution dne to nnlt angles of attaok extending Inward from the tip . 
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lift and drag components so as to maintain them along 
the perpendicular and the parallel to the actual local 
relative-wind velocity. A further refinement of the 
theory would involve the influence of the curvature of 
the wing wake. Since the helix angles involved in the 
motions are small (p6/2V<0.1 and r&/2F<0.1) and 
since that region of the wake nearest the wing is of 
predominant influence, this correction may be neglected. 

Inasmuch as the various stability derivatives thus 
depend upon a summation of appropriate components 
of the lift and the drag loading along the span, it was 
necessary to determine these distributions for each of 
the wings with several different angle-of-attack distri- 
butions. For this purpose the Lotz method of calcula- 
tion (references 1 and 2) was used. In order to keep 
the computations from becoming too lengthy, the chord- 
distribution function that occurs in this method was 
expressed by, at most, three terms of a cosine series (as 
in table I). Although this expression caused the chord 
distributions of the actual wings (fig. 1) to differ slightly 
from those for linearly tapered wings with rounded tips, 
such a procedure was justified because these slight de- 
partures in plan form had only a small effect on the 
characteristics but permitted a large saving in the com- 
putations required. Thus only the terms near the 
diagonal running through equations (19) of reference 2 
entered into the computations. As the various deriva- 
tives for the elliptical wing could be obtained relatively 
easily, they were sometimes computed in order to deter- 
mine the shape of the various derivative curves; it was 
therefore possible to use fewer points in fairing similar 
curves for the tapered wings. 

The wing theory was applied in a special way so as 
to obtain results applicable to any arbitrary twist of 
the wings. The theoretical span loading being a linear 
function of the angle-of-attack distribution, the loading 
due to arbitrary twist can be built up, as will be indi- 
cated later, from certain elementary loadings by super- 
position. The elementary loadings considered were 
those caused by simple unit jumps of angle of attack 
occurring at different points of the span. 

For each of the nine tapered wings, the first 20 
Fourier coefficients doterrnining the load distributions 
were computed (10 odd and 10 even) for the cases of 
unit angles of attack extending inward from the wing 
tip and covering various amounts of the semispan. 
The rest of the wing was in each case assumed to be at 
zero angle of attack. The portions thus covered were 
0.25, 0.50, 0.75, and 1.00 of the semispan. 

In spite of the great number of harmonic terms re- 
tained, the conditions near the points of discontinuity 
in the angle of attack required special treatment. The 
problem of these end conditions has been solved by 
Betz and Petersohn (reference 3) and their results were 
utilized in fairing the load curves through this region. 
Figure 2 shows the elementary loadings that were cal- 
culated, including the modified fairing. The results 
pertain specifically to the chord distributions illustrated 



in figure 1 but interpolation might be made for inter- 
mediate plan forms. 

It is evident that any angle-of-attack distribution, 
symmetrical or unsynimetrical, may be built up of 
elemental steps of the type used in deriving figure 2. 
Figure 3 illustrates the procedure of finding the resultant 
load distribution. Thus, the loading contributed by 
element- 10 of figure 3 (a) is obtained by deducting tho 
load curve due to an increment of angle of attack extend- 
ing between B and C from that due to an increment 
extending between A and C. Although this process 
could be continued until the load distribution was com- 
pletely determined, the same results can be more easily 
obtained from influence lines, which give tho load at a 
particular spanwise station due to tho effect of unit 
angle-of-attack changes extending inward various 



(a) 




A-C. 



(a) Angk-of-attaok distribution. 
Cb) Load distribution Tor element (10). 
Floras 3.— Load components lor an olemont. 

amounts from the light wing tip. Such influence lines 
are given in figures 4, 5, and 6 for eight evenly spaced 
points across the wing semispan. Each lino was ob- 
tained by cross-plotting the values of at tho inter- 
sections of the loading curves of figure 2 with vertical 
lines drawn at the particular stations. (For example, 
the points of intersection of line Z — Z, fig. 2 (b), with 
the various curves represent the load induced at the 
0.25-semispan point by uniform angle-of-attack incre- 
ments that extend in varying amounts from the wing 
tip. These intersections identify the corresponding 
curves of figs. 4 to 6.) 

In order to illustrate tho use of the influence lines in 
determining the lift distribution as well as to show the 
degree of accuracy with which they may be U6ed, the 
influence lines will be applied to predict tho loading for 
a tapered wing (X=0.25, A=&) corresponding to the 
angle-of-attack distribution shown in figure 7 (a). 
The particular angle-of-attack distribution used is 
defined by the equation 



(c, 



- . i ■ " ) sin 30 
4o sin 8/ 



(3) 
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Relative distance from wing center 
(a) X.-L00. 00 X-0.S0. C<0 X-0J5. 

Ficcbe 4.- Influence lines for determining the load distribution, X=5. The nnmber Identifying a given line refers to the particular spanwise station at which the load 

is to be computed. 
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-1-0 -3 -.6 -.4 -.B O .2 A .6 .8 I.Q 

Relative distance from wing center 



(a) X-1.00 (b) X=0JS0. (o) X=0.2S. 

5,— Influence lines for determining tbe load distribution. A- 10. Tire number Identifying a given line refers to the particular spanwlse station at which the 

load Is to be computed. 
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Relative distance from whg center 
(a) \=IM. (b) X=0J0. Co) X=0.25. 

Fiol-re 8— Influence lines for detetminini the load distribution. A-=16. The number identifying a given line refers to the particular spanwlse station at which the load 

[3 to be computed. 
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This particular distribution is employed because it is 
possible thereby to compute exactly the corresponding 
theoretical distribution as a check, without the usual 
approximations of a Fourier series. The following pro- 
cedure illustrates the use of the influence lines to deter- 
mine the lift at the 0.75-semispan point due to this dis- 
tribution of twist: (1) The influence curve labeled 
0.75 in figure 4 (c) is reproduced beneath the angle dis- 
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Fiouee 7.— Determination of lift distribution. 



1.0 



tribution to the same spanwise scale; (2) a base line 
with a range from a max - to a ma3!+ is laid out as in figure 
7 (c) with the origin of the qrdinates at a equal to zero; 
(3) the effect of any length of elemental angle-of-attack 
change, da, in figure 7 (a) is found by projecting the 
length of the element onto figure 7 (b) and plotting the 
increments (A{) and (A 2 +A 2 ') at the angles of attack 



for which these elements are drawn, as in figure 7 (c)- 
Because a negative angle would induce a negative load 
at the point in question, Ai is plotted as a negative 
value. This process is continued from a MM _ to a MX |. 
and the resulting curve (fig. 7 (c)) is integrated to obtain 
the total effect at 0.75, which is then plotted in figure 
7 (d). The load distribution over the entire span is 
obtained by repeating the same procedure for a number 
of points along the span. 

With the lift loading thus determined, the induced- 
drag distribution may be found by a simple operation, 
namely 



(4) 



Figure 7 (d) gives the comparison of the load- 
distribution curve obtained from tho influence lines 
with that computed directly by the wing theory using 
equation (3). Although the agreement is not precise, 
it must be remembered that tho solid curve represents 
a case where no series approximation was necessary; 
hence it may be concluded that tho influenco-lino 
method of determining the lift distribution is as ac- 
curate as any other for practical purposes. 

Aside from other possible applications, tho load dis- 
tribution may be used to determine the stability deriva- 
tives for certain cases not specifically covered by the cal- 
culations. In the subsequent charts, it is sometimes 
necessary to stipulate either that the initial anglo-of- 
attack distribution be symmetrical about the wing 
center line or that the dihedral angle be constant along 
the span. With a knowledge of the complete load 
distribution, however, values of the derivatives or their 
respective moments might be found for particular cases 
where the charts do not apply. 

STABILITY DERIVATIVES 

Although it is possible, in the general case, to obtain 
the stability derivatives from the lift distribution, such 
a procedure will not usually be necessary because the 
charts to be presented cover all cases likely to be of 
interest. The results are presented in such a form that 
the effect of flaps on the derivatives may easily bo 
determined. 

ROLLING MOMENT DUE TO ROLLING 

The first derivative considered is the rolling moment 
due to rolling. In unstalled flight, when the wing rolls 
about the longitudinal wind axis, a damping or restoring 
moment is set up. This moment L„ lliH , varies directly 
with the angular velocity p and is defined by the 
equation 



(5) 



where the product C^^pis simply a rolling-moment 
coefficient that varies linearly with the angular velocity. 
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The computed variation of the derivative Ci p 'with 
aspect ratio and taper is given in figure 8. In the usual 
lift range below the stall, this derivative may be con- 
sidered to be independent of either initial wing twist or 
angle of attack and of the wing dihedral. In conven- 
tional cases, practically the entire damping moment for 
the airplane may be attributed to the wings. It can be 
seen that the moment contributed by a tail surface, 
geometrically similar to the wing but with only one- 
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Aspect ratio, A 
FiGCBI 8.— Rolling derivative dua to rolling. Z„m»i =Ci, |p.jS6 

fourth the span, would be Km of that contributed by 
the wing, inasmuch as L TOl i ins for a given angular ve- 
locity varies as the fourth power of a linear dimension. 

Reducing the aspect ratio or increasing the taper 
tends to reduce the derivative C^, as may be seen from 
the curves given in figure 8. Comparison of the pres- 
ent values of with similarly derived values given in 
reference 4 indicates that the effect of rounding the 
wing tips is to reduce the theoretical restoring moment 
by about 6 percent for wings of aspect ratio 6. 

YAWING MOMENT DUE TO ROLLING 

During a rolling motion, the wing experiences a linear 
antisymmetrical change in angle of attack along the 
span and, as a result, antisymmetrical loadings are 
added to those that originally were on the wing. The 
resulting yawing moment is due to components of the 
lift as well as to the drag along the span, the lift com- 



ponents being the more important. With the specified 
system of axes, positive rolling produces a negative 
yawing moment or, for any case, with positive lift coef- 
ficients the falling wing tends to advance owing to the 
predominating influence of the lift vectors. 

The yawing moment due to rolling, unlike the rolling 
moment due to rolling, depends upon both the initial 
wing twist and the angle of attack. For untwisted 
wings, however, the yawing moment is zero at zero lift 
and increases linearly with the wing lift coefficient. 
For a twisted wing, the yawing moment due to rolling, 
although varying linearly with the over-all lift coeffi- 
cient, is not necessarily zero when C L is zero but may 
have either a small positive or a small negative value 
depending upon the initial angle-of-attack distribution. 
Owing to tliia circumstance, it is most convenient to 
express the derivative C ap as a ratio in terms of unit 
partial-span angle-of-attack changes. 

Figure 9 shows the computed variation of the ratio 
Cnja. for unit symmetrical angle-of-attack changes that 
extend out from the wing center so as to cover various 
amounts of the wing span. Thus, if it is desired to 
determine C ap for an untwisted rectangular wing of 
aspect ratio 6 at an angle of attack of 0.1 radian, the 
value 0.195 (for an angle of attack of 1 radian), read 
from the solid line of figure 9 (a) at the relative dis- 
tance of 1.0, is multiplied by 0.1 to give a value of 
C„ equal to 0.0195. If, now, a half-span flap of con- 
stant chord ratio were displaced an amount sufficient 
to- cause an additional change in angle equal to 0.1 
radian over the portion with flaps, the new value of 
would be 

(0.1X0.195)+(0.1X0.134)=0.0329 1 
This value of is then inserted into the equation 

CL^ffSJ (6) 

to determine the yawing moment due to a rolling 
angular velocity. 

Although the curves given in figure 9 can be directly 
used to determine the effect on G Kp of deflecting partial- 
span flaps of constant flap-chord ratio, they are also 
readily adapted to the determination of C Kp for a wing 
with any initial twist provided that the twist distribu- 
tion is symmetrical about the wing center line. The 
process is illustrated in the following example where it 
is desired to find the value of 0 Kp for a rounded-tip 
rectangular wing of aspect ratio 6 with the symmetrical 
angle-of-attack distribution shown in figure 10 (a). 
The contribution of the element of angle of attack da, 
shown at the point a=0.15 radian, to the total value 
of the wing C Sp is equivalent to that caused by a full- 
span elemental flap min us the contribution of the 
cross-hatched portions. The contribution of this 

element da is denoted byf— Vee and may be obtained 
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(a) A -6. 

<b) A =10. 
(o) A-U. 

Figuhe 9.— Yawing derivative dae to rolling. 

by projecting the various small flap portions onto the 
appropriate C n Ja curve (taken from fig. 9) as in figure 



10 (b), and adding the increments A! and A a . The 
sum of these increments is then plotted in figure 10 (c) 
at the value of a for which the element is drawn. The 
value of C Up for the complete wing is obtained by per- 
forming the integration 



(7) 



8 
t. 



yOTm 


1 














1 




















il 
























eft* 








i 
i 






n 












1- 

1 
I 






i 

i 

— i. 
i 
i 
i 












T 
1 
1 







Refofii/a distance 



C np 













r 

i 
l 






r — n 
1 
i 












h 

l 
l 

1- 




-a 


r 

i 

>)4 
i 
i 












i 
i 
i 












i 










































l 


h 






Root 


















Tip 



.6 



.8 



from wing center 



1.0 



O .2 .4 .6 .8 

Extent of unif symmetrical angle of attack 



IX) 









































>) — 


















— (« 




































\ 








Cn, ' 


• 0.0 


391 


* 




















( 




5 



.05 JO .15 .20 

ct. radian 

Figure 10.— Application of C, f curves to an oxamplo. 



.25 



These curves apply to wings with symmetrical twist 
and it is necessary to consider only half the wing, the 
factor 2 being included in the curves. The evaluation 
in the case of figure 10 (c) yields 0.0391. 
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The curves of figure 9 indicate that, for a given 
angle-of-attack distribution, there is relatively little 
change in the value of C„ v 'with the taper ratios investi- 
gated. Changes in taper ratio did, however, have an 
appreciable effect on the value of (fig. 8) with the 
result that the ratio of the yawing to the rolling moment 
in roll will, in general, increase with increase in taper. 

Inasmuch as the inclination of the lift vectors at the 
outer portions of the span has such a predominating 
effect on the yawing moment, the most effective means 
of reducing for a given wing lift coefficient is to 
give the wings washout toward the tips. 

The yawing moment due to rolling is, in conventional 
designs, largely due to the wings. The tail surfaces 
contribute very little to this moment both because of 



approximations in deriving the necessary equations for 
the determination of the yawing derivatives. TVhen 
these approximations are used and the velocity along 
the span is expressed as a variable, the new downwash 
equation becomes 

w =±r P- dB v 2T „ } . 2T ke (8) 

4t6 Jo , cos 8— cos 6 0 

The system of simultaneous equations derived for the 
approximate solution of this integral equation is 

ZCa, cos 2nd 2 A. sin nB +^2nA K sin nB 

iP 2 (^-t-^+i) sm n6 = ~ 5 » si* n9 = a si n 0 (9) 
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Figure II.— Rolling derivative doe to yawing. 



their short span and because of the small angles of 
attack relative to the wing. 

ROLLING MOMENT DUE TO YAWING 

During a yawing motion, increments of velocity are 
added along the forward-moving half of the wing and 
similar increments are deducted along the rearward- 
moving half. The difference in velocity of the two 
halves causes a rolling moment which, for an untwisted 
wing, varies directly with the initial angle of attack as 
well as with the angular velocity. The velocity incre- 
ments vary linearly with the distance from the wing 
center line and are small relative to the flight speed; 
it is therefore permissible to make certain mathematical 



in contrast to the system given by equation (18) of 
reference 2. 

By means of equation (9), Fourier coefficients were 
computed for the nine tapered wings with two different 
initial angle-of-attack distributions: (1) a distribution 
due to a unit angle of attack extending over the whole 
span, and (2) a unit angle of attack at the wing center 
covering half the span. In order to obtain the correct 
fairing of the final curves of figure 11, similar results 
were computed for elliptical wings with six angle-of- 
attack distributions covering 0, % % %, %, and all of 
the wing span. 

As was the case with the derivative C Bp , it is most 
convenient to give the derivative of rolling moment 
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due to yawing O tr as a ratio in terms of a partial-span 
unit angle of attack. The values of C tr may be ob- 
tained from figure 11 and are to be inserted into the 
equation 

rb 

Lya«ta,=Gi t ^yqSb (10) 

By the process described in the previous section, values 
of C tr may be obtained for wings with any initial twist 
distribution that is symmetrical about the center line. 



.6 




O .2 .4 .6 .8 W 

Extent of flap from iv/ng center 



(a) A =6. 

(b) .4=10. 

(c) A=16. 

Figure 12.— Yawing derivative duo to yawing for partial-span flap. 

The curves of figure 11 fall in the order that would 
be expected for the various taper ratios, i. e., the 
moment for an untwisted tapered wing would be ex- 
pected to be less than that for a rectangular wing of the 
same span and area because the tapered wings have a 
smaller proportion of the wing area at the tip. On 
account of the induced velocities along the span, the 
reduction, for the tapered wings, is not so great as 
would be obtained by an application of the ordinary 
strip theory. 

The direction of the moment is such that, with the 
system of axes used, a positive rolling moment generally 



results from a positive yawing velocity when the wing 
is giving positive lift. By the use of considerable 
washout, such as is obtained with partial-span flaps, 
it is possible not only to reduce the value of this moment 
but also to make it slightly negative for low wing lift 
coefficients. 

As was the case with and C*,, the value of C, r 
for the entire airplane is duo almost wholly to the 
wings because the side area of the airplane contributes 
relatively little moment as compared with the wings 
in curvilinear flight. 
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Floras 13.— Yawing derivative dua to yawing tor untwisted wing. 



TAWING MOMENT DUE TO TAWING 

A part of the wing yawing moment due to yawing 
results from the change in the induced-drag distribution 
that accompanies the change in the lift distribution 
across the span of a yawing wing. The rest of the 
yawing moment is due to the difference in the distribu- 
tion of profile drag resulting from the variation in 
velocity along the span. Both parts, however, pro- 
duce damping moments in the unstalled-flight range. 

The part of the wing yawing moment duo to the 
induced drag is defined by the equation 

N veiata ,=C AT ^.q.Sb (11) 

wherathe derivative C, r may be obtained from figure 12 
for certain types of angle-of-attack distribution. Tho 
special distributions for which the derivatives of figure 
12 apply are both uniform and symmetrical about the 
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wing center line. Such distributions occur only when 
partial-span flaps of constant-chord ratio are deflected, 
the rest of the span being at zero angle of attack. This 
limitation in the applicability of these curves as com- 
pared with the previous ones is due to the fact that the 
principle of superposition does not apply in cases where 
the variation is not linear with a. The computed 
results may, however, be used to determine the varia- 
tion of C n , for the most useful case, namely, that of a 
wing without twist. For this purpose, the proper 
values of C K J<£ obtained from figure 13, which is a cross 
plot of the end points of figure 12, are multiplied by the 
square of the actual angle of attack. 

The part of the yawing moment due to the profile 
drag can be determined from the easily derived equation 

6 

&N M =zf\ cj(l +^}f dy (12) 

where c*,, and c are functions of the distance y along the 
span. It is possible, by assuming constant and by 
neglecting terms of the second order, to obtain a coef- 
ficient A& r that may be used with the equation 

*N vav{ *,=AC« T ^qSb (13) 

to compute the part of the yawing moment due to the 
profile drag. The values of the profile yawing mo- 
ment, as given by equation (13), are sufficiently 
accurate for most wings since Cs a generally varies only 
slightly across the span. The variation of the coef- 
ficient ACxJCdq with taper ratio is given in figure 14. 

•4\ 1 1 1 1 1 , , , , , 




0 .2 .4 .6 .8 /.O 

Taper rath, A 

Figure It— Proflle-drag yawing derivative. &N,.* l „-i.C^^iSt 

The total wing yawing moment due to yawing is the 
sum of the moments given by equations (11) and (13). 
At low lift coefficients, the profile drag contributes the 
greater portion of the wing damping moment in yawing. 
At moderate or high lift coefficients, however, the part 



due to the induced drag exceeds that due to the profile 
drag. If it is assumed that a=0.3 and (7^=0.01, then 
the respective values of G ar and AC Xr would be 0.0522 
and 0.0031 for a rectangular wing of aspect ratio 6. 

The damping moment contributed by the wings in 
yawing motion is, in most cases, secondary but is not 



Leading edge---. 




Resulting angle-of-aftack distribution 
Fiousz 15— Effect of Irregular dihedral on sideslip. 



negligible with respect to the damping moment con- 
tributed by the fuselage and the vertical tail surfaces. 
The damping in yawing due to the wings depends upon 
the angle of attack as well as upon the plan form; 
therefore the relative amounts contributed by the 
wings and tail surfaces may vary considerably. 

Although it was not possible to give a general chart 
for detennining the damping in yawing for symmetri- 
cally twisted wings as was done with the previous deriva- 
tives, it can nevertheless be said that the addition of 
load toward the tips, whether by washin or by an 
increase in taper ratio, would increase the wing damping 
moment due to a yawing angular velocity. 

ROLLING MOMENT DUE TO SIDESLIP 

The manner in which the changes in angle of attack 
that cause a rolling moment are brought about during 
a sideslipping motion is shown in figure 15 by a sketch 
of a wing having positive, negative, and zero dihedral 
over various portions of the span. For simplicity, the 
wing is assumed to have no initial twist and the dihedral 
angles are assumed constant over each of the portions 
A, B, and C. For small angles of sideslip /3, the in- 
crease in angle at tip A is, to a first approximation, 
equal to r^j3; whereas, at the opposite tip A', there is 
an equal decrease of the angle of attack. The portions 
B — B', having no dihedral, contribute no change in angle 
of attack when the wing is sideslipping. At the center, 
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however, owing to the negative angle of dihedral T c , 
there is an effective decrease in angle of attack over 
part C equal to TcP and on C there is a similar increase 
in angle. Figure 15 shows the resulting effective 
angle-of-attack distribution for the particular shape of 
dihedral assumed. 

The effect of this distribution is similar to that caused 
by two pairs of ailerons equally and oppositely deflected 
with the inner pair opposing the rolling action of 
those at the tip. Positive areas of dihedral on the 
advancing wing tend to add load onto that wing. For 
the system of axes chosen, all areas with positive 
dihedral produce a negative rolling moment with a 
positive angle of sideslip. This moment, like the rolling 
moment due to roll, is independent of the initial wing 
twist as long as no portion of the wing becomes stalled. 

The rolling-moment derivative due to sideslip 
may be determined from figure 16, which gives the 
variation of CijY for various unit antisymmetrical 
angle-of-attack distributions (i. e., symmetrical portions 
with constant dihedral) that extend out from the wing 
center and cover various relative amounts of the wing 
semispan. In the usual case, where the dihedral 
angle T is constant along each semispan, the value of 
the rolling moment due to a sideslip angle /3 can be 
obtained from the equation 

OyftSb (14) 

where the appropriate values of C^fT, obtained from 
figure 16 at the relative distance equal to 1.0, are multi- 
plied by the dihedral angle in radians. In more 
unusual cases as, for example, where only the tips are 
turned up or where the wing is given a gull shape for 
any reason, it is still possible to determine a coefficient 
of rolling moment due to sideslip simply by adding the 
effects of the various parts in the way previously de- 
scribed. Thus, for the wing shown in figure 15, let 
A=6, X=1.0, T A and r c =0.1 radian and assume that 
it is desired to find the proper value of (7^ to use in 
equation (14). The part due to the tip portions A — A' 
is A, (from fig. 16 (a)) X IU=0.195 X 0.1=0.0195. 
The part due to the center portions C — C is A 2 X T c = 
0.065X0.1=0.0065. The resulting value of C, fi to 
be used in equation (14) is thus 0.0130. The extension 
of this method to a curvilinear variation of r along 
the span may be easily made by plotting the values of 
T at each point of the span and using the method given 
in a previous section for integrating for the total effect. 

The results of figure 16 indicate that equivalent 
angle-of-attack changes caused by unit lengths of 
dihedral portion near six-tenths of the relative distance 
from the center are, in general, slightly more effective 
in producing rolling moment than . unit lengths of di- 
hedral at the tips. Such a result is due partly to the 
fact that the load curves near the tips are rounded and 
partly to the fact that, for the tapered wings, the larger 
areas affected by lengths of dihedral near the 0.6 point 
tend to compensate for the shorter moment arms 
through which the change in loading acts. 
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Although, during a sideslipping motion, positive 
dihedral produces a righting moment, a similar though 
generally smaller effect may also be produced by the 
addition of vertical area above the longitudinal wind 
axis. Also, on account of interference effects, the 
proportion of the airplane rolling moment contributed 
by the wings may vary considerably with the external 
appearance of the airplane. 

It is usually considered, in practice, that a straight 
wing will have some dihedral effect, but tests of wings 
with well-rounded tips (reference 5) do not support this 
view. In cases of wings with blunt tips or in cases 
where chords of the sections near the tip do not lie in 
one plane, some dihedral action is shown. 
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FiGiai 17— Comparison between experimental and computed values of CiffT (ex- 
perimental data tram reference 5). 

Figure 17 shows a comparison of experimental and 
computed values of C^/r. The experimental values 
have been obtained from figure 23 of reference 5 and the 
coefficients given therein have been converted to the 
form used in this report. In the tests reported in 
reference 5, a rounded-tip rectangular wing of aspect 
ratio 6 was given various lengths of dihedral by turning 
up the outer portions of the wing. Each wing was then 
tested throughout the angle-of-attack range for various 
sideslip and dihedral angles. 



It will have been apparent from the preceding dis- 
cussion that the results of figure 16 may also be applied 
to predict the rolling moment caused by an aileron 
deflection in unyawed flight since ailerons, equally and 
oppositely deflected, cause changes in the angle-of- 
attack distribution that are similar to the changes 
caused by dihedral. Strictly speaking, however, the 
change in angle of attack due to dihedral cannot have 
quite the same effect as a similar change produced by 
ailerons because the ordinary lifting-line theory, when 
applied to yawed or sweptback wings, omits the effect 
of the stagger of the trailing vortices and the inclination 
of the bound vortex. Although the present theory has 
not been modified to take this effect into account, there 
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Flop or- aileron chord/ Wing chord 
Figcbx 18.— Variation of k with ratio of flap or aileron chord to wing chord. 



is ample justification for omitting it in the computa- 
tions as experiments indicate only second-order differ- 
ences (see reference 5) for the usual angles of yaw and 
sweepback. 

For the computation of the rolling moment due to an 
aileron deflection 5, the appropriate value of C lt to be 
inserted in the equation 

Ci&Sb (15) 

may also be found from figure 16. The derivative C fJ 
is given as a ratio in terms of k, the theoretical change 
of a with aileron deflection. Although the value of k 
has been theoretically determined for thin wings, it is 
better to use values of k determined from an analysis 
of experimental data. For this purpose, figure 18 is 
included, which shows the variation of k for values of 
the ratio of aileron or Sap chord to wing chord up to 
0.3. This variation of k has previously been given in 
figure 1 1 of reference 6 and holds for sealed flaps deflected 
up to approximately 20°. 
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(a) Ailerons deflected equally. 

(b) Ailerons deflected differentially. 
Figure IS— Addition of effect of aileron elements. 
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c, qcT/3dy 




Front view 

Figure 20.— Action of dihedral In producing yawing moment In sideslip. 

If the angle-of-attack change caused by deflecting 
the ailerons is antisymmetrical about the wing center 
line, the proper value of O t> to be used with equation 
(15) (for the rolling moment only) can be found by an 



integration or summation of the effects of elemental 
ailerons of various lengths and positions along tho span 
as indicated in figure 19 (a). The values of C tl {k are 
obtained from figure 16 for the wing plan form used. 
If the ailerons are differentially operated, then it may 
be better to divide the ordinates of figure 16 by 2 and 
to determine the value of the moment given by each 
aileron as indicated in figure 19 (b). 

YAWING MOMENT DUE TO SIDESLIP 

The yawing moment of a wing with dihedral in side- 
slipping motion may be conveniently divided into two 
parts, the first part being due to the unsymmetrical 
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FiotraK 21.— Yawing derivative dne to sideslip (dihedral constant). 
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induced-drag distribution over the span and the second 
part due to a shift of the lift vectors acting so as to pro- 
duce a moment about the vertical axis. Figure 20 illus- 
trates the components of the section lift and drag 
vectors that produce yawing moments. The advanced 
wing having the larger lift will also have a larger induced 
drag and hence a component moment is set up that tends 
to turn the wing so as to reduce the sideslip ; at the same 
time, however, the contrary moment due to the com- 
ponents of the lift acts to advance the forward half of 
the wing still more. As was the case with the yawing 
derivative due to rolling, the moments caused by the lift 
components predominate and, as a result, the net 
theoretical moment is an unstable one; or, in other 
words, with the system of axes chosen, a negative 
yawing moment results when the dihedral and sideslip 
angles are positive. 

The explanations advanced in some textbooks neglect 
the inward slope of the lift vectors and lead to an incor- 
rect sign of the yawing moment. 

The yawing moment in sideslip is given by the equa- 
tion 

N. i *e, li ,=C lt0 0< l Sb (16) 

The derivative C np is given in figure 21 as a ratio in 

terms of Ta because its value depends linearly upon the 
magnitude of the product of these variables. The 
values of 0„ JV a have been computed for unit symmetri- 
cal angle-of-attack distributions that extend out on 
either side of the center line and cover 0.25, 0.50, 0.75, 
and all of the wing span. These carves may be used to 
determine values of C K for any initial angle-of-attack 
distribution symmetrical about the wing center line, 
provided also that the angle of dihedral is constant 
across the wing span. Although the rolling derivative 
due to sideslip can be obtained (from fig. 16) for a curvi- 
linear variation of dihedral along the span, it is necessary 
to stipulate that either a or r remain constant if the 
principle of superposition is to be applied in the deter- 
mination of C„ g . The combination of variable symmet- 
rical twist and uniform dihedral being more co mm on 
than the converse, the computations were shortened by 
including curves for only the case of uniform dihedral. 

The resultant value of C np (to be used in equation 
(16)) is found by either an integration or a s umm ation 
of the effects of elements of angle of attack extending 
along the span. The process to be followed where 
graphical evaluation is necessary has been illustrated 
in figure 10, with the ordinates of figure 10 (a) changed 
to Ta. The ordinates and abscissas of the remaining 
parts are to be changed as required. For untwisted 
wings with uniform dihedral, the value of C^Ta. is 
obtained by multiplying the value read at a relative 
distance of 1.0 by the wing angle of attack and, in 
turn, by the dihedral angle. 

The curves of figure 21 being generally steeper be- 
yond the 0.5 point, the deduction of increments of 



angle of attack at the tip, i. e., giving the wing wash- 
out, would be the simplest means of decreasing the 
unstable yawing moment caused by the wings in a 
sideslipping motion. 

Although the predicted variation of the yawing 
moment with dihedral is confirmed, experiments show 
a residual stable yawing moment at zero dihedral that 
is not predicted by the ordinary theory. This residual 
moment is greater for wings with blunt tips and is 
greater at zero or negative lifts. It will be noted that 
the theoretical yawing moment is itself the small re- 
sultant of two large contrary effects and is thus of the 
same order as a number of possible secondary influences. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., April 19, 19S8. 
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Table L— VALUES OF COEFFICIENTS DEFINING WING 
CHORD DISTRIBUTION 
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